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Within the Relativistic Theory of Gravitation it is shown that the equation of state 
p = p holds near the center of a black hole. For the stiff equation of state p = p — p c 
the interior metric is solved exactly. It is matched with the Schwarzschild metric, which 
is deformed in a narrow range beyond the horizon. The solution is regular everywhere, 
' with a specific shape at the origin. The gravitational redshift at the horizon remains 

C ~ ) ' finite but is large, z ~ 1O 23 M0 /M. Time keeps its standard role also in the interior. The 

QQ energy of the Schwarzschild metric, shown to be minus infinity in the General Theory of 

f***^ , Relativity, is regularized in this setup, resulting in E = Mc 2 . 

Keywords: Black hole interior, stiff equation of state, vacuum equation of state, Rela- 
tivistic Theory of Gravitation, bimetric theory, exact solution 



1. Introduction 

C3 



Black holes (BH's) have fascinated mankind because of their predicted properties: 
at the classical level not even light can escape; the mass is located in the center; 
the role of time in the interior is played by space and vice versa; on the quantum 
level they radiate as a thermal body via the Hawking mechanism. 

The structure of galaxies is related to astrophysical BH's, galaxies are believed 
to have a supermassive central black hole with mass equal to 0.12% of the bulge 
mass. [1] The BH in the center of our own Galaxy has "only" four million solar 
masses. Some believe that our Universe is actually the inside of a giant black hole. 
[2] 

As the above mentioned theoretical aspects are physically difficult to understand, 
an important question is: Are they perhaps unrealistic, and should the behavior 
of physical BH's not be sought in a modified version of the General Theory of 
Relativity (GTR)? One such candidate is the Relativistic Theory of Gravitation 
(RTG) proposed by Logunov. [3,4] Following idea's of Rosen [5] and others, e.g. 
[6,7], it writes the Hilbert-Einstein Lagrangian in terms of a field in Minkowski 
space-time, extends it with the cosmological term — p\ = —8irGA/c 2 and with a 
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bimetric coupling \pbiliiv9^ v between the Minkowski metric 7 and the Riemann 
metric g. [3, 8] The bimetric term breaks the general coordinate invariance of 
GTR, and allows only the harmonic gauge for the metric or the Lorentz gauge 
for the field. Since it will be cosmologically small, [4] RTG will have the same 
content as GTR for all standard observable effects in the solar system and galactic 
problems. Differences may arise near singularities, in particular, for black holes and 
in cosmology. 

An important question is the value and, in particular, the sign of the new pa- 
rameter phi- Logunov pointed out that the choice pbi = PA cancels the zero point 
energies of both terms, allowing to have far away from matter just a Minkowski 
space. When RTG was formulated, good data for p\ were not available, so it was 
natural to choose /?bi < 0, with a graviton mass m g = — p bi /f 671. To make up 
for the observed positive value of the cosmological constant, an inflaton field can 
be added. [4] In [8] we have considered the situation where pbi = PA > 0, its value 
being set by the present cosmological data, making an inflaton field obsolete. This 
leads to a tachyonic graviton, but its tachyonic nature sets in only at the Hubble 
scale, where not individual gravitons but the whole Universe is relevant. 

The bimetric coupling regularizes the infinite redshift at the horizon of black 
holes [3,4]. In a previous approach we presented a scaling behavior near the horizon: 
coming from the outside, the time-time and radial-radial components of the metric 
tensor follow the Schwarzschild shape, but they cross over to an exponential decay 
in the interior. [8] They do not change sign, and thus leave for time its standard 
role, which is physically appealing: BH's are then extreme objects, but still behave 
similar to normal ones. Indeed, with the Killing vector remaining time-like in the 
interior, the Hawking mechanism does not work, as it neither does for Newton stars. 

In [8] the solution in the interior was considered on the basis of present consensus: 
all matter in or very near the origin. We modeled this by a very-low-pressure 
equation of state. In our follow up studies, we have realized that this approach is 
inconsistent. Here we report about the opposite case: matter spreads throughout 
the BH. It is modeled by the stiff equation of state p — p — p c . In GTR this shape 
is known to have simplifying features [10] and this appears to carry over to RTG. 

2. Setup of the problem 

Static spherically symmetric bodies have a metric 

ds 2 = U(r)c 2 dt 2 - V(r)dr 2 - W 2 (r)(d6 2 + sin 2 ^ 2 ). (1) 
and lead in RTG to the 00 and 11 Einstein equations 

1 W' 2 2W" VW _ 8irG 
W 2 ~ VW 2 ~ VW + V 2 W ~ c 4 PtoU 
1 , W' 2 U'W _ 8ttG 
W2 + yW2 + JJ^F ~ — Pt°t' ( ) 

where the total density and pressure have the form 

. . Pbi Phi PhiT 2 

. Pbi Pbi . PbiT 2 , - 

Ptot = P-P^ + ^j-2V + -^- ( 3 ) 
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The harmonic constraint reads 



U' V W' 4rV 

h 4 = — T . (4) 

U V w w 2 y ' 



Taking G = c = 1 we define the mass function M (r) by 

W' 2 

V= 1-2M/W (5) 
The Einstein equations can now be written as 

We shall neglect p\ and, because U <C min(F, 1), only keep the pbi/U terms. 
Indeed, they are responsible for the scaling behavior near the horizon. [3, 8] The 
fundamental assumption of the present Letter is that they also determine the shape 
in the interior of the BH. 

To start, let us neglect matter and suppose that U ~ W 2 . Eq. ([6]) first brings 
A4 ~ W and then sets U — 8Trpt,iW 2 . The energy conservation condition reads 

(p + P )U' + 2p'U = 0. (7) 

If we assume that p and p are bounded at small W, then the result p' ~ —W'/W 
leads to a logarithmic divergence, in conflict with the assumed boundedness. The 
case p = np implies p ~ W^ 1 ^ 1 ^. For every k (0 < n < 1) this is more singular 
than the presumed leading term pbi/ (2f ) ~ W~ 2 . Thus for a singular RTG solution 
one cannot treat matter as a perturbation, rather one must take p ~ p at large p, 
that is, near the origin. 

Let us therefore consider the stiff equation of state [9] 

p = p-p c , r < R; p = p = Q, r > R, (8) 

From energy conservation there arises the shape 

p = l Pc {Jk-i), p = i Pc (^ + i), (9 ) 

where the subscript c denotes the point r = R. 

Notice that the vacuum equation of state, p = —p is covered in case they are 
homogeneous. It appears in the limit U c — ► 0, x c — > 0. 

3. Exact solution for the interior 

The problem now leads to a similar exact solution. In terms of 



W / 3 

1 - iv Wl = Vss- (10) 



the result J7 = 87r(p c [/ c + ph\)W 2 may be written as 



U = U c X - 2 , x c =J P f c v (ID 

£c V SlPc^c + Pbi 
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while density and pressure read 



/< = ^c(?| + l), P=^>c(|t-1)- (12) 



We may define at the point x = 1 

E! 

P (X: I) 

This allows to express 



p c (l-3af|) 2(1 + 2ki) 

Both factors in the left hand side being positive (p c because p < p; U c to avoid a 
horizon {/ = 0), it is seen that the sign of pbi is set by the physically allowed value 
of K\. For "classical" matter, it is natural to assume that negative pressures do not 
occur, and that p — at the horizon, so k.\ — 0. This is the case to be considered 
from now on. (Quantum matter may have k\ down to —1; We come back to this in 
the discussion.) The bimetric and cosmological coupling constants are the assumed 
to be negative. We now have 

The harmonic constraint thus brings 

x' x" 2xx' x' 8rx' 2 

2 ^- 2 V-—* + ^ = xHi^-y (16) 

Going to the inverse function r{x) makes it linear, 

x 2 (l- x 2 )r" + x(3- Ax 2 )r' = 4r. (17) 
Let us define the conjugate variable 



y = Vl-x 2 , (18) 
The solution is then remarkably simple, 



r = n (l + ^=)x^ 1 (l + y)-^ 7 (19) 

where r% is the value at x — 1. Near that point one has r = ri(l — -^y)- We can 
now derive from (fH)]) and (fTO|) . 



4r*i 

V = ^ x 4-2V5 {1 + y) 2V5, (2Q) 
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Our solution is thus completely explicit. At the origin it exhibits the singularities 
known for the stiff equation of state, [4, 10] 

U = U x r^, V = -flWfr't*- 2 , W = W x r^», (21) 

where 7 M = |(v5 + 1) is the golden mean. But if we take W as the coordinate, we 
have the Ricmann metric 

ds 2 = _^£w2 dt 2 2d ^ _ WW, (22) 

W x 2 x% 1-W 2 /4:M 2 ' v ; 

in the interior of the BH. It clearly is regular at its origin, with the factor 2 coding 
the above singularities, and a coordinate infinity (but not a change in signs) at the 
horizon. 

4. Matching with the exterior 

Well away from matter, the harmonic constraint brings for W the Schwarzschild 
shape Ws = r + M, where M = Ai(R) is the mass, basically equal to the mass as 
observed at infinity. This implies via Eq. (T5|) 

Vs = l - = L±^i = ±. (23 ) 

1 - 2M/Ws r-M U s 

In GTR one has pbi — 0, so Eq. gives x c = 1/V3. Taking W as coordinate 
sets W — ► 1 in (J5J), yielding V c = 3 = 1/U C , W c — 3M. This solution is seen as a 
limit model for neutron stars. [10] But let us consider the matching problem in the 
present setup at some r = R, x = X . Since we have to match W' = 1, it follows that 
V = 2/(1 -X 2 ), implying X 2 = \{R+ZM)/{R + M). On the other hand, equating 
M = i(l + X 2 )Ws to M, yields as only solution R = M, X = 1, exceeding the 
presumed maximum 1/V3- Thus, in the harmonic gauge the set (V, W, W) cannot 
match the Schwarzschild values, so this limit model of GTR must be distrusted. 

In general, we consider as regular any solution for which M(R) < -kW(R). We 
define a black hole as a solution for which M(R) « hW(R). For a black hole Eq. 
(fT5)l sets x c = 1 and Eq. (fTl"]) shows that this is possible within RTG, provided pbi 
is negative, which is the Logunov situation with a massive graviton. Together with 
Eqs. fTO]). |H]) this amounts to 

W > = 2M ' r ' = a (24) 

Contrary the Schwarzschild philosophy, we demand that V remains finite, as it 
occurs in (|38[) at y — 0. Eq. ([5]) then offers an equivalent manner to characterize a 
BH, 

Criterion for black hole horizon: W (R) ~ 0. (25) 

It is handy to introduce the inverse length m g ("graviton mass") and the di- 
mensionless small parameter rn g , 
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where we took the estimate from Ref. [4] . We shall therefore have the values at the 
horizon R 

K Jlf2 

U c = m% V c =- 1 p, W C = 2M, W' C = Q. (27) 

These values look worrying, as they are far from Schwarzschild's, even when r is 
near M (e.g., W' s = 1). The problem nevertheless appears to be consistent. Beyond 
R we need the deformation of the Schwarzschild metric near the horizon, caused 
by the bimetric coupling. It was studied by Logunov and coworkers [3,4], and an 
elegant scaling form for small m g was presented by us, [8] 

r = M ] + f+ll r °\ , U = ^, (28) 
1 -ry(e« + £ + r ) 

e« 2M 
V = — ^, W = 



r/(l + e«) 2 ' 1 -7?(e« +w ) +m 2 g £' 

Here £ is the running variable, rj a small scale and r$ and wq parameters. Coming 
from the outside, = 0(l/rj), the functions follow the Schwarzschild shapes, with 
small corrections, but they branch off for = 0(1). At £ = the function V 
has a maximum l/(4ry), while U has already gone down to r\. Going further to the 
inside, for £ <C —1, U and V both decay exponentially over a very short distance, 
Sr = 2r\M . It is this exponential decay that will provide the opportunity to match 
the seemingly very different behaviors near the horizon. Matching U, V and W 
with the boundary values (|27[) of the interior solution, we find that 

wo = \[\ w h fiogrnfl + \ log h - - lV (29) 



In this regime Eq. (f2"B"]) yields 



W\r) = ±— f — -f . (30) 



7y(e« + 1) 

The values (|2"9"|) confirm that W'(R) = 0, at the considered order m g . This property 
nicely settles a subtlety. The scaling shape of W becomes negative below a certain 
r, [8] a fact erroneously interpreted as self-repulsion. [3,4] If matter is taken into 
account, it induces a W > in the interior, which goes to zero at the horizon. This 
matches the zero coming from the outside. All by all, one thus has W > 0, except 
for W'(R) = 0, see Fig. 1. 

To fix the parameters ro and wo of Eq. (|28|) , we need to see which effects they 
bring at finite r. Around the Schwarzschild solution (neglecting for now ph\ and 
Pa) there are four perturbative modes. The first, 

(rL-M)M rL~M LM 2 

°Ui = — 7 — , ,„-,„ , oVi - 



2 ' 



(r + M) 2 ' r-M (r - M) 

1 1 r + M 

SWi = -(ri — M), L = -iog-Z_, (31) 
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involves the logarithm of Vs — l/Us- The second mode, (5U2, SV2, 5W%) = d(Us,Vs, 
Ws)/dM, relates to a shift in the mass; the third, (£/s,0,0), rescales Us, while the 
fourth, (0, Vg, hWs), rescales Vs and Wg. We require that the second is absent at 
order m g and, to keep the proper behavior at infinity, that the third and fourth 
are absent to all orders. At order m g ~ 77 this may be imposed by analyzing the 
behaviors of U, V and W of Eqs. (|28[) for r > M, with r — M small but fixed, 
yielding 

r = 2 + logri, w = O{m g ), (32) 

the latter being in agreement with (|29|) . The logarithmic mode (|3Tj) remains with 
prefactor 877, and it is coded in the terms linear in £ of Eq. ([28]) . It is a finite 
distance, nonperturbative effect of order m g ~ M^—pbi- Near the horizon, for 
r—M ~ fn g M ln(l/m s ) , it signifies the onset of the deformation of the Schwarzschild 
metric. 




Fig 11. The metric function W(r) has value 2M and slope zero at the horizon of a black hole. 
m = 0.01. 



" log(U) / r/R 
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Fig 12. The metric function U starts as a powerlaw and becomes equal to nig at the horizon, be- 
yond which it grows exponentially towards the Schwarzschild shape. m g = 1.5 10 -23 corresponds 
to a one solar mass black hole. 
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Fig 13. The metric function V is of order unity in the interior, with a powerlaw divergence at the 
origin. Beyond the horizon it grows exponentially towards a maximum of order l/m g , after which 
it joins the Schwarzschild shape. m g = 10 -3 . 



From Eqs. (J2HJ) and (f2U)) we get the horizon radius, 



R = M 



l + A\j-m g (log m g + 1) 



< M. 



(33) 



With the metric completely specified, we present plots of U, V and W in Figs. 
1, 2, 3, respectively. 

Having determined the leading scaling approach, we may go to next order in r\ 
in the peak regime = 0(1) and in the horizon regime where it is 0(m g ). This 
rather painful analysis will not be reported here. We just mention the confirmation 
of our black hole condition W'(R) — at second order in m g . 

All by all, we may now rewrite the scaling form by eliminating £ in favor of £/, 



= M 



V 

w 



l + [/ + m gV /275(log£/ + 2) 
m g y275(logC/ + 2)' 



U - 
U 



(34) 



1 - U + m 2 + mi \og(U/m 2 ) 



This describes the free space region r > R, where rn 2 g < U < 1 + 0(m g ). At scale 
r ~ 1/m Newton's law is picks up the Yukawa factor exp(— mr), due to the massive 
nature of gravitation in RTG. [3] 

The interior shape can also be expressed in U as running variable, where it lies 
in the range (0,m^). Due to Eqs. (JTTJ) and (TT5J) it also holds that 



U 



V = 




(35) 
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The density and pressure read 



1 ttl 2 1 in 2 

P= 2 Pc{ -U +1) > P =2 Pc{ -U- 1} > (36) 



With R given by Eq. (O, the locus is 



r = R(l + - yf^^Hl + y)-^ +1 \ (37) 

v5 

and the other two metric functions read 

V= 5 ^(l-yf-^(l + yr +V ~ 5 , W = 2M^. (38) 

The behaviors of the metric functions U, V and W are plotted in Figs. 1, 2 and 
3, respectively. 

5. Properties of the solution 

First of all, with p and p' , the functions U', V , W" are discontinuous at horizon. 
For U' this is possible in the Eq. (2), because W'(R) = 0. 

The characteristic size of the deformation range of the Schwarzschild solution, 
^deform = r(£ = 1) - R, is small, 



e 1+e \/2 M 2 

Reform < \l ^m g M In « 1.6 HT 18 — j m. (39) 

V 5 m„V5 Ml 



For one solar mass BH's this is comparable to the Compton radius of the W and Z 
bosons, i w = 2.45 lCr 18 m. 

The gravitational energy density was discussed elsewhere. [8,11] For the metric 
((!]) it takes the form 



i00 c A W 2 ( r 2 V'WW 3t . r 2t . 

8nGr 6 \ V 



2r 3 VW 
W 



8rWW' - 2r 2 V - 3W 2 . (40) 



At the origin it diverges as 5 , which is integrable. H The total energy density 
reads 6 00 = t m + VW 4 ptot/r' 1 - [8] Its separate contributions are depicted in Fig. 
2. 



In static, spherically symmetric RTG the gravitational energy density is unique, "gravitational 
energy can be localized" . This is implied by the harmonic constraint. The residual gauge transfor- 
mations allow only shapes that are unbounded at or oo. As they would make the energy infinite, 
they have to be discarded. This holds both for regular solutions (stars) and for our BH solution. 
2 A different energy momentum tensor was proposed in Ref. [7]. Its energy density diverges at 
the origin as r(^~ 9 )/ 2 , too singular to give a finite integral. For this reason we shall abandon it. 
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We can calculate the material energy. A partial integration is needed to nu- 
merically tame the divergent behavior near r = 0, that is exposed in Fig. 4. This 
brings 

f R o VW 4 

U mat = 4tt / drr 2 — —p = 2228.830945 M, (41) 
Jo r 

which is pretty large. Gravitational terms are negative and subtract 2227.830945 M 
from this. The bimetric term in p t ot contributes as 

?7bi = 47r [ drr 2 ^ = -1372.93286 M. (42) 
Jo r 2(7 

The gravitational energy inside the BH is 

f R 

(7 graVjint =/ dr47rr 2 i 00 = -842.898079 M. (43) 
To 

Together they make up for (7i n t rior = 13 M. 

The gravitational energy density in the skin layer first has a large positive and 
then a large negative part, due to the term r 3 V , see Fig. 4 in the region around 
r/R = 1. The integrated effect is obtained easily since the formulation of the 
Einstein equations in Minkowski space implies that the total energy density is a 
total derivative, 



1 d /VW 2 + W 4 W 3 W 
Airr 2 dr \ 2r 2r 3 



0°° = 7^^( 1 ^- + ^T- 1 -T-)- (44) 



where we have set c = G = 1. In the Schwarzschild regimeff] so its this combines 
into 

B°° = _L ± M{r + Mf{2r + M) 

Airr 2 dr 2r 3 (r — M) ' { ° ! 

From flU}, (051 and (271), tne region i? < r < oo yields [/exterior = M(l- 5-8 + 0) = 
— 12M. Together with the interior it makes up the total BH energy U = Mc 2 . 

6. Conclusion 

We have considered a black hole in the Relativistic Theory of Gravitation. Previous 
findings that the bimetric coupling regulates the divergencies of the Schwarzschild 
singularity, are extended to show that it sets the behavior in the interior. Near 
the center, consistency requires that p fts p is unbounded. For the case of the stiff 
equation of state an exact and rather elegant solution is provided for the interior. 
It matches the deformed Schwarzschild solution of the exterior. For a BH of one 
solar mass, the deformation range is of the order of the Compton length of the W 
and Z bosons. 



3 For the Schwarzschild problem (matter-free Einstein equations in GTR), in the standard gauge 
W s = r, U s = l/Vs = 1 - 2M/r, we obtain 00 = -M 2 /{2nr 2 (r - 2M) 2 ]. Like 05), it has a 
quadratic divergence at the horizon, which is non-integrable and poses a so far overlooked problem 
for the Schwarzschild metric of GTR. 
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Fig 14. Full curve: The total energy density © 00 as function of r/R for m g = 0.00375. With 
V' , it is discontinuous at the horizon. For realistic situations, very small rn 9 , the peak near the 
horizon is much higher, narrower and deeper on its right side. Dashed: Material energy density 
VW 4 p/r 4 . Dots: Gravitational density t 00 + V W 4 p^i/2Ur 4 in the interior. Dash-dots: Energy 
density of the Schwarzschild metric; in RTG its divergence at the horizon is regularized by the full 
curve. 



Powerlaw singularities occur at the origin. This has been a reason to discard 
the problem, [4] but they disappear when the standard radial coordinate W(r) is 
employed, rather than r itself. Away from the origin, the solution is regular, and in 
particular also at the horizon, for any observer. The redshift at the horizon is finite, 
though of the order l/rn g > 10 23 M Q /M. In the interior, time keeps it standard 
role. Hawking radiation is absent, and Bekenstein-Hawking entropy has no bearing. 

Open problems are to derive the radiation and to treat, for a given type of 
matter, the equation of state self-consistently with the metric, as was done here near 
the origin for any type of matter. Next, quantization of the field theoretic approach 
can be considered. One may also extend the approach to the Kerr-Newman black 
hole. 

Though we have elegantly described the interior of a black hole, we have not 
been able to settle definitively the question of the sign of the bimetric coupling 
Pbi- Indeed, while Eq. (fl4|) definitely leads to pbi < for non- negative pressures 
(«i > 0), we cannot yet exclude the regime —1 < m < —h, where a positive pbi 
would be required. Let us mention in this connection that the situation with the 
vacuum equation of state p — —p {n\ = — 1) holding in the interior was recently 
analyzed and stated to describe a "gravastar" in its Bose-Einstein condensed ground 
state. [12] We plan to investigate the connection between RTG black holes and Bose- 
Einstein condensation in the near future. In that situation values of k\ near — 1 
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may occur, in principle. Let us mention that our exact metric remains valid for the 
vacuum equation of state, that is, in the limit Ki — > — 1, U c — > 0, and x c — > 0. 

Having a complete solution at hand, we could also verify some general aspects 
of the gravitational energy momentum tensor. We realized that the energy density 
of the standard Schwarzschild metric of GTR has at the horizon a quadratic sin- 
gularity. The related infinite gravitational energy is not apparent in the Riemann 
approach, the singularity then being viewed as a coordinate singularity. [13] When 
making the step from GTR to RTG, the divergence at the horizon gets regular- 
ized [8], and here we have seen that the total energy is finite, and equal to Mc , 
as it should. We also verified that for static, spherically symmetric bodies, residual 
gauge transformations within this sector would lead to infinite energies. So they are 
forbidden, making the local energy density uniquely defined: energy and, in par- 
ticular, gravitational energy, can be localized in RTG, probably under more general 
conditions than reported here. 

The fact our black hole has no true horizon would justify as name: grey hole. 
However, if observed black holes in the cosmos have a huge but finite redshift at 
their horizon, as described here, it is better to stick to the standard name. 

The resolution of the singular behavior at the horizon arises from the bimetric 
coupling, which acts as a mass-type term. This breaks the general coordinate 
invariance of GTR. Since pbi is cosmologically small, this could normally play a role 
only in cosmology - but it still allows the A Cold Dark Matter model. Indeed, it 
brings no change of general relativistic effects in the solar system or for gravitational 
radiation of binaries. However, we have seen that the bimetric term does play a 
role at large redshifts, that is to say, near the horizon of Schwarzschild black holes 
and inside it. This resolution of a singularity may be more general. 

Returning to the black hole problem: R is sometimes argued that our Uni- 
verse may actually be the inside of a giant black hole. [2] For that application, the 
Schwarzschild metric with all its matter in the center is not realistic, and our setup 
with matter spread throughout the interior looks more natural. 
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